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' Abstract 



We consider the random walks killed at the boundary of the quarter plane, with 
homogeneous non-zero jump probabilities to the eight nearest neighbors and drift zero 
in the interior, and which admit a positive harmonic polynomial of degree three. For 
these processes, we find the asymptotic of the Green functions along all infinite paths 
Q\ | of states, and from this we deduce that the Martin compactification is the one-point 

compactification. 
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+2 ! 1 Introduction and main results 

c3 

First introduced for Brownian motion by R. Martin in 1941, the concept of Martin 
compactification has then been extended for countable discrete time Markov chains by J. 
Doob and G. Hunt at the end of the fifties. The purpose of this theory is to describe the 
asymptotic behavior of the Markov chains and also to characterize all their non-negative 
I/" - ) \ superharmonic and harmonic functions, see e.g. [Dyn69 1 . 

For a transient Markov chain with state space E, the Martin compactification of E is 
the smallest compactification E of E for which the Martin kernels z *— > k z z ° = G z ° /G z z l 
extend continuously - by G z z ° we mean the Green functions of the process, i.e. the mean 
| number of visits made by the process at z starting at zq, and we note z\ a reference 

state. E \ E is usually called the full Martin boundary. Clearly, for a E E, zq h-> k§> is 
superharmonic ; then d m E = {a € E \ E : zq i— > /c^° is minimal harmonic} is called the 
minimal Martin boundary - a harmonic function h is said minimal if < h < h with h 
harmonic implies h = ch for some constant c. Then, every superharmonic (resp. harmonic) 



function h can be written as h(zo) = k z ° fi(dz) (resp. h(zo) = J 9 E k z fi(dz)), where \i 
is some finite measure, uniquely characterized in the second case above. 

The case of the homogeneous random walks in is now completely understood. First, 
their minimal Martin boundary is found in [DSW60], thanks to Choquet-Deny theory. 
Furthermore, in the case of a non-zero drift, P. Ney and F. Spitzer find, in their well- 
known paper [NS66], the asymptotic of the Green functions, by using exponential changes 
of measure and the local limit theorem ; this gives consequently the concrete realization 
of the Martin compactification, in that case the sphere. Additionally, in the case of a drift 
zero, the asymptotic of the Green functions is computed in |Spi64| ; it follows that the 
Martin compactification consists in the one-point compactification. 

Results on Martin boundary for non-homogeneous random walks are scarcer and more 
recent. We concentrate here our analysis on important and recently extensively studied 
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examples that are the random walks in X, d killed at the boundary of cones. They are 
related to many areas of probability, as e.g. to non-colliding random walks or quantum 
processes. 

On the one hand, the case of the non-zero drift is now rather well studied. 

In |Bia92b] . P. Biane considers quantum random walks on the dual of compact Lie 
groups and, by restriction, arrives at classical random walks with non-zero drift killed at 
the boundary of the Weyl chamber of the dual. Solving an equation of Choquet-Deny 
type, he finds the minimal Martin boundary of these processes. 

When the compact Lie group is SU(d) and the associated random walk has non-zero 
drift, the Martin compactification is obtained in ColO 1[. by finding the asymptotic of the 
Green kernels. 

Recently, in [I R09] . I. Ignatiouk-Robert obtains the Martin compactification of the 
random walks in with non-zero drift and killed at the boundary. She uses there an 
original approach based on large deviations theory in order to compute the asymptotic of 
the Martin kernels. Unfortunately, her methods seem quite difficult to extend up to the 
case of the drift zero. Also, they do not provide the asymptotic of the Green functions. 

This asymptotic in the case of the dimension d = 2 is found in [KR09J . 

On the other hand, results on Martin boundary for killed random walks with drift 
zero are quite rare. The simplest example of the cartesian product is due to [PW92] . A 
more interesting case comes again from quantum processes : in [Bia92aJ, P. Biane shows 
that the minimal Martin boundary of the random walk with zero drift and killed at the 
boundary of the Weyl chamber of the dual of SU(d) is reduced to one point. 

It is immediate from |Bia91j that this classical random walk in the Weyl chamber of 
the dual of SU(d) is, for d = 3, the random walk spatially homogeneous on the lattice 
{i+j exp(z7r/3), € Z 2 } with jump probabilities as represented on the left of Picture[TJ 




Figure 1: Random walk in the Weyl chamber of the dual of SU(3) 

Obviously, a suitable linear transformation maps the lattice {i + j exp(z7r/3) : € 
Z 2 } into Z 2 , see Pictured]; in this way, the Weyl chamber {i + j exp(z7r/3) : € 
becomes Z+. For d = 3, the killed random walk considered by P. Biane in [Bia91j can 
therefore be viewed as an element of 

= {random walks in Z+ with non-zero jump probabilities ipij)—i<ij<i 

to the eight nearest neighbors, with drift zero and killed at the boundary} 

with jump probabilities as represented on the right of PictureQ]- above, "drift zero" means 
that pi,i+pi,o+Pi ,-i = o+p-i - i andp^i+po.i+P-i.i = Pi ,-i+Po -i- 

In this setting, P. Biane proves, in [Bia92a], that (io, Jo) ^ Wo^o+io) is the only positive 
harmonic function for this process. 

By the same methods, it can certainly be shown that there is only one positive 
harmonic function for the "dual" walk, namely for the random walk with jump probabilities 
P-l,o = Po,i = Pi,-i = 1/3. In particular, if we set & c = {random walks of such that 
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Po-l = P-1,1 = Pi,o = V, P-i,o = Po,i = Pi ,-l = v , V + v = V 3 } " in otli er words, & c is 
the set of all cartesian products of the random walk on the dual of SU(3) with its dual, 
see on the left of Picture [2] below — , it follows from [PW92] that any process of £P C has 
also a minimal Martin boundary reduced to one point. 
In this paper, we introduce the set 

=^1,0 = {random walks of for which (io, Jo) ^ *ojo(*o + jo) is harmonic}. 

Note that we have £P C C &i,o, but we will see, in Remark HI that the inclusion is strict. 
More generally, we define 

&a,p = {random walks of for which (io, Jo) ^ Wo(*o + Q Jo + 0) is harmonic}. (1) 

Since any harmonic function for a killed process takes the value zero on the boundary, 
&a,p is in fact exactly the set of all killed random walks in to the eight nearest 
neighbors admitting a harmonic polynomial of degree three. 

The description of the set & a ,p i n terms of the (pij)i,j is rather cumbersome but not 
difficult to obtain, it is postponed until Remark [U Let us just note here that if a > 2 or 
a < 1/2, then for all /3, = ; if a = 1/2 or a = 2, then for all /3 / 0, = 0, 

and & a ,o is reduced to one walk ; and if a G]l/2, 2[ and is small enough, then ^ a ,p 
is a (non-empty) set with two free parameters, properly described in Remark [U We have 
represented on the right of Picture [2] an example of a process belonging to & a ,0i f° r an Y 
a €[1/2,2]. 

Moreover, note that considering in this paper t^afi is all the more natural as the set 
{random walks of for which (io,Jo) l— ^ *oio is harmonic} is studied in [Ras09j. 
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Figure 2: On the left, a generic walk of 2? c {fx + u = 1/3) ; on the right, an example of walk 
of ^ Qi0 (A = a(a-l/2)/[2-a+2a 2 ], // = (a/2)/[2-a+2a 2 ] and v = (l-a/2)/[2-a+2a 2 ]) 

Our first result deals with the Green functions - below, (X, Y) denotes the coordinates 
of the random walk and E(j j ) the conditional expectation given (X(0), Y(0)) = (io, jo) - 



.fc>0 



J (*oJo) 



{(X(k),Y(k))=(i,j)} 



(2) 



and, more precisely, with their asymptotic along all paths of states. 



Theorem 1. Suppose that the process belongs to & a ,f3- Then the Green functions (0) 
admit the following asymptotic as i + j — > oo and j/i — >• tan(7), 7 lying in [0, 7r/2] : 



Gtf - Ciojo (io + a Jo + 0) W^l> , (3) 



+ aj) 
[i 2 + mj + a 2 j 2 

where C > depends only on the parameters (pi,j)i,j and is made explicit in the proof 
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In the particular case of the random walk killed at the boundary of the Weyl chamber 
of the dual of SU(3), the asymptotic §3§ is, for 7 G]0,7r/2[, proved in |Bia91j . Theorem[T] 
actually completes this result for that very particular random walk and, in fact, gives the 
asymptotic of the Green functions for a much larger class of processes. 

In addition, Theorem Q] has the following consequence. 

Corollary 2. The Martin compactification of any process belonging to &a,/3 is the one- 
point compactification. 

Furthermore, the asymptotic ([3j) of the Green functions in the two limit cases 7 = 
and 7 = ir/2 enables us to obtain the asymptotic of the absorption probabilities 

K' 30 = V {io>jo) [3k>l:(X(k),Y(k)) = (i,0)], 

hf j ° = r (ioJo) [3k>l:(X(k),Y(k)) = (0,j)]. W 

Indeed, the absorption probabilities are related to the Green functions ([2]) through 

h i ~ V\ + P0 

hj — P-l.lCrj j_ 1 +P-l t0 U 1 j +P-1 -lCr lj - +1 , 

so that, from Theorem [H we immediately obtain the following result. 

Corollary 3. Suppose that the process belongs to 8? a fi. Then the absorption probabilities 
0) admit the following asymptotic as i — > 00 : 

K M ~ C(Pi,-l +P0-1 +P-i,-i)«oio (io + ujo + 0) -tj, 

where C > is the same constant as is the statement of Theorem^ 

The same asymptotic holds for h ®' 10 , after having replaced (pi _i+po — — l) above 
by (p-1,1 +P-i,o +y-i,-i)/a 5 . 

The asymptotic of the absorption probabilities in the case of a non-zero drift being 
obtained in [KR09J, Corollary [3]thus gives an example of the behavior of these probabilities 
in the case of a drift zero. 

In order to prove Theorem [lj we are going to develop methods initiated in [FIM99j and 
based on complex analysis, what will allow us to express explicitly the Green functions ([2]). 
Indeed, in [FIM99] . the authors G. Fayolle, R. Iasnogorodski and V. Malyshev elaborate 
a profound and ingenious analytic approach for studying the stationary probabilities for 
random walks to the eight nearest neighbors in the quarter plane supposed ergodic, i.e. 
such that 271,1+231,0 +P1-1 < P-1,1 +P-i,o +P-1 ,-1 and Pi,i +Po,l + P-1,1 < Pl,-1 +Po,-l + 
P-1,-1- 

We are going to see here that this analytical approach can be extended up to the case 
of the random walks in the quarter plane with drift zero and killed at the boundary : 
Section [2] of this paper first broadens the analysis begun in Part 6 of [FIM99] for the drift 
zero, and then shows how this applies in the case of the random walks of £P a a. 

It is worth noting that this approach via complex analysis is intrinsic to the dimension 
d = 2 ; for this reason, it seems really a difficult task to generalize it in higher dimension. 

Let us conclude this introductory part by describing the set defined in ([T|) in 

terms of the jump probabilities (pij)i,j- 

Remark 4. The fact that the two drifts are equal to zero gives two equations and the fact 
that the sum of the jump probabilities is one yields an other one. Moreover, the harmonicity 
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of h(i ,j ) = i j (i + aj + /3), which reads h(i ,j ) = ^ijPi,jh(i + i,jo + j), leads to 
ten equations, by identification of the coefficients of the third-degree polynomials above. It 
turns out that some of these equations are trivial and that some other ones are linearly 
dependent, we finally obtain six equations linearly independent. We can therefore express 
all the eight jump probabilities (pij)i,j in terms of p\,\ and pi^ only, and we obtain : 

* P-1,0 = -["(1 ~2a- 13) + 8pi,i + (4 - 3a + 2a 2 + a/3) Plfi ]/[a{l + 2a + /?)], 

* p_i fl = [a(l - a - p) + 2(4 + 3a + 2a 2 + a/3) Pl<1 + 2(2 + a 2 + a(3)p lfi ]/[2a(l + 2a + 0)}, 

* P0,l = -[-(1 + a + P) + 4(2 + 2a + P) PlA + 2(2 + a + /3)pi, ]/[2(l + 2a + $)], 

* pi _i = [a 2 + (-1 + 2a - - (1 + /3 + 2a 2 )pi, ]/[l + 2a + /3], 

* Po-l = -[(-l-3a-/3+4a 2 )+4(-2+2a-/3)p l5l + (-4+6a-2/3-8a 2 )pi i0 ]/[2(l+2a+/3)], 

* p_!_a = [a(l-3a-/3+2a 2 )+2(4-3a+2a 2 -a/3)pi i i+2(2-3a+3a 2 -2a 3 )pi j0 ]/[2a(l + 
2a + '/3)]. 

The properties of ^a :/ 3 mentioned below ([1]) are immediately obtained by studying the 
sign of the jump probabilities above in terms of a, /3, p% t i and pi^. 

2 Explicit expression of the Green functions 

Section [2] aims at obtaining an explicit expression of the Green functions ([2]) - what we will 
succeed in doing in Theorem [8] below. This forthcoming expression of the Green functions 
will be, in turn, the starting point of Section [3j where we will find their asymptotic. 

In order to prove Theorem El we need to state two results, namely Equation ([6]) and 
Proposition [6] : Equation (|6|) is a functional equation between the generating function of 
the Green functions (|2|) and the ones of the absorption probabilities ([!]), and Proposition [6] 
establishes some quite important properties of the generating functions of the absorption 
probabilities. 

The proof of Proposition[6]turns out to require considering the Riemann surface defined 
by {(x,y) € (CU{oo}) 2 : Yli j Pi,j x% D^ = 1}> f° r this reason we begin Section[2]by studying 
- and, in fact, by uniformizing - this surface. 

It seems of interest to us to introduce this Riemann surface in whole generality ; this 
is why, at the beginning of Section [21 we are going to suppose that the process belongs to 
0* - and not necessarily to S?a,$- 

To begin with, we define the generating functions of the Green functions ([2]) and of 
the absorption probabilities (jlj) by 

G io ' jo {x,y) = Gj/V-y- 1 , h io ' jo (x) =Y J h i ?' jo x\ h io ' jo (y) = ^// : - , V (5) 

and 4°o° = F (t ,j )i 3k ^ 1 : ( x ( k ), Y ( k )) = (0,0)]. Of course, G io ^°, h io ^° and h io 'i° are 
holomorphic in their unit disc. With these notations, we can state the following functional 
equation on {(x,y) € C 2 : \x\ < 1, \y\ < 1} : 

Q (x, y) G io ^° (x, y) = h io 'i° (x) + h io ^° (y) + hftj - x io y>° , (6) 

where Q(x, y) = xy [ ^ a Pi,jX l y : ' — l] . Equation ([6]) is obtained exactly as in Subsection 2.1 
of [KR09| . 

The polynomial Q(x,y) defined above can obviously be written as 

Q(x, y) = a(x)y 2 + b(x)y + c(x) = a(y)x 2 + b(y)x + c(y), 

with 

a(x) = pi,ix 2 + p ,ix + p-i,i, b(x) =Pi,qx 2 -x + p-i fi , c(x) = p\ -\x 2 +p -ix + P-\ _i, 
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a(y) =pi,iy 2 +pi,ox +pi -i, b(y) =p ,iy 2 - y + po-i, c(y) = P _ 1A y 2 + p-i fi y + p-i-i- 
Let us also define the polynomials 

d{x) = b{x) 2 - 4a(x)c(x), d(y) = b{yf - Aa{y)c{y). 

It is proved in Part 2.3 of [FIM99] that for any random walk of d (resp. d) has one 
simple root in ] — 1, 1[, that we call x\ (resp. yi), a double root at 1, and a simple root in 
E U {oo} \ [—1, 1], that we note X4 (resp. 2/4). 

For example, in the case of SU(3), i.e. for the random walk with transition probabilities 
as in Picture [TJ we immediately obtain x\ = 0, y% = 1/4, x A = 4 and 7/4 = 00. 

From a general point of view, it is shown in Part 2.3 of [FIM99| that x\ (resp. yi) 
is positive, zero or negative depending on whether p-ifl 2 — 4p_i j ip_i ) _i (resp. po,-i 2 — 
4pi ! _ip_i j _i) is positive, zero or negative, and that X4 (resp. 7/4) is positive, infinite or 
negative depending on whether pi : o 2 — ^Pi : iPi,-i (resp. Po,i 2 ~~ 4p 11 p_ 11 ) is positive, zero 
or negative. 

Let us now have a look to the surface defined by {(x,y) € (C U {oo}) 2 : Q(x,y) = 0}, 
that we note for the sake of briefness. Note first that Q(x,y) = is equivalent to 
[6(x) + 2a(x)y] 2 = d{x) or to [b(y) + 2d(y)x] 2 = d(y). As a consequence, it follows from the 
particular form of d or of d (two distinct simple roots different from 1 and one double root 
at 1) that the surface has genus zero, and is thus homeomorphic to a sphere C U {00}, 
see e.g. Parts 4.17 and 5.12 of [JS87]. Therefore, this Riemann surface can be rationally 
uniformized, in the sense that it is possible to find two rational functions x(z) and y(z), 
such that i2 = {(x(z),y(z)) : z € C U {00}} ; moreover, a standard uniformization (for an 
account of the concept of uniformization, see Part 4.9 of [JS87] ) is : 

_ ( z - Zl )( z -l/ Zl ) _ (z-Kz 3 )(z-K/z 3 ) 

1 ] ~ (z-z )(z-i/ Zo y y[z) ~ (z-Kz 2 )(z-K/ Z2 y {n 

where 

z = [2-(x 1 +X4) + 2[(l-x 1 )(l-x 4 )] 1/2 ]/[x4-x 1 ], 

Z\ = \xi + £4 — 2xiX4 + 2[xiX4(l — Xi)(l — X4)] 1 / 2 ] / [X4 — Xi] , 

Z2 = [2-(y l +y A ) + 2[(l-y l )(l-y A )} 1 ' 2 ]/[y A -y l ], 

Z3 = [yi +y± - 2yiy A + 2[yiy 4 (l - y x )(l - y 4 )) l/2 ]/[y A - yi], 

and where K is a complex number of modulus 1. Note that zq and z\ (resp. z 2 and z 3 ) have 
a modulus equal to one or are real, according to the signs of xi and X4 (resp. y± and y A ). 
For example, in the case of SU(3), it follows from a direct calculation that 

zq = exp(— 2ot/3), z\ = 1, z 2 = exp(— 27r/3), 23 = exp(«7r/3), K = exp(— z7r/3). 

Above and throughout the paper, we note 1 the usual complex number verifying i 2 = — 1. 

In the general case, in order to find K, we need to introduce a group of automorphisms 
naturally associated with the surface 2. To begin with, let us remark that, with the 
previous notations, Q(x,y) = entails Q(x,[c(x)/a(x)]/y) = and Q([c(y)/d(y)]/x, y) = 
; it is therefore natural to consider the group generated by the two bilinear transfor- 
mations £(x,y) = (x,[c(x)/a(x)]/y) and fj(x,y) = ([c(y)/a(y)]/x, y), which is called, in 
[FIM99] , the group of the random walk. 

These automorphisms f; and f) define two automorphisms £ and n of the uniformization 
space C U {00}, characterized by : 

£o£ = 1, xo£ = x, yo£ = [c(x)/a(x)]/y, 7707? = 1, yorj = y, xor] = [c{y)/a{y)]/x. (8) 
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With ([7J) and flSJ), we obtain that they are equal to : 

= 1/z, V (z) = K 2 /z. (9) 

In particular, it is immediate that the group W = (£,ri) generated by £ and 77 is isomorphic 
to the dihedral group of order 2inf{n > : K 2n = 1}. For example, in the case of SU(3) 
for which K = exp(— iir/3), W is of order six - this fact is (differently) proved in Part 4.1 
of |FIM99j . 

A crucial fact is that this property is actually verified by any random walk of & a R, 
since we have the following. 

Proposition 5. For any process of ^ a ^, K = exp(— m/3). 

Proof. With (0), we have y(K) = y l ; in addition, by ©, rj(K) = K, so that with ©, 
we obtain x(K) 2 = c(yi)/a(y{). This implies that x(K) = — [c(yi)/a(yi)] 1//2 - indeed, we 
easily show that the roots of Q(x, y\) have to be negative. By using again ([7|), we get 
K + l/K = [a{ yi ) l l 2 { Zl + 1/zi) + c(y 1 ) 1 /2( ZQ + i/ Zq )] / [ a { yi ) l l 2 + c{ yi ) l l 2 ] . In particular, 
K + 1/K can be expressed explicitly in terms of the jump probabilities (pij)i,j- By using 
then Remark [Hand after simplification, we get K = exp(— nr/3). □ 

From now on, we suppose that the process belongs to & a ,p- 

For a better understanding of the surface £} as well as for a coming use, we are now 
going to be interested in the transformations through the uniformization (x,y) of some 
important cycles, namely the branch cuts [xi,X4], [2/1,2/4] and the unit circles {\x\ = 1}, 
{|y| = 1}. First, by using (|7|) and Proposition [5j we immediately obtain : 

z- 1 ([xi,x 4 ]) =MU{oo}, y- 1 ([yi, 2 / 4 ]) = exp(-z7r/3)RU{oo}. (10) 

As for the cycles rc _1 ({|x| = 1}) and 2/ -1 ({|j/| = 1}), their explicit expression (calcu- 
lated starting from ([7])) shows that they are real elliptic curves, which are located as in 
the middle of Picture [3] below. 



1*1=1 




Figure 3: The uniformization space C U {00}, with on the left some important elements 
of it, in the middle the corresponding elements through the uniformization (x,y), and on 
the right the images of the cone F = {xex.p(i9) : x > 0, — 7r/3 < 6 < 0} through the six 
elements of the group W = (£, rf) 

Note also that with Q and Proposition [SJ we immediately obtain £(exp(2#)lR + ) = 
exp(— z#)lL|_ and n(exp(i9)R + ) = exp(—z(9 + 27r/3))M + . In particular, if we denote by F 
the set {xexp(i6) : x > 0, — ir/3 < 9 < 0}, we have - see also on the right of Picture [3]- 

(jw(F) = C. (11) 
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Thanks to the group W = (£,77) and to (fTT|) . we are now going to continue the lifted 
functions H io >i°(z) = h io ' jo (x(z)) and H io ^(z) = h io ' jo (y(z)) ; this fact will turn out to 
be of the highest importance in the proof of Theorem [8]- the latter being crucial, since it 
will be the starting point of the forthcoming Section [3j 

Note that in the sequel, we are often going to write x l °y 3o {z) instead of x{z) l °y{z) 30 . 

Proposition 6. The functions H io >i°(z) = h io > jo (x(z)) and H io ^°(z) = h io > jo (y(z)) can be 
meromorphically continued from respectively {z G C : \x(z)\ < 1} and {z G C : \y(z)\ < 1} 
up to respectively C \ exp(z-7r)[0, 00] and C \ exp(2«7r/3)[0, 00]. These continuations verify 

H io,jo ( z ) = gio jo > H iojo = H io,jo ( v ( 12 ) 
for all z E C, and 

if z E C\D (13a) 

■ (-1)' H x io y jo («>(*)) if zeD (13b) 



H io ' jo (z)+H io ' jo (z)+h^°-x io y jo (z) 



where we have set D = {xexp(iO) : x > 0, 2tt/3 < 8 < it} and l(w) for the length of w, 
i.e. the smallest r for which we can write w = w\ - ■ ■ w r , with w\, . . . ,w r equal to £ or rj. 

Remark 7. In {z G C : \x(z)\ < 1, \y(z)\ <l}cC\fl, \13a\) follows directly from (0|). 

Proof of Proposition In order to prove Proposition [6l we are going to use strongly the 
decomposition (jlip : precisely, we are going to define H 10 ' 30 and _£P 0j0 piecewise, by 
defining them on each of the six domains w(F) that appear in the decomposition (jlip . 
to be equal to some functions Hw and Hw' 30 ■ It will then be enough to show that the 
functions H Mh3 ° and H 10,30 so defined verify the conclusions of Proposition [6) 

• In F = {xexp(i9) : x > 0, -vr/3 < 9 < 0} C {z G C : \x(z)\ < 1, \y(z)\ < 1}, see 
Picture [H we are going to use the most natural way to define H l0 ' J0 and H 10,30 , i.e. their 
power series. So we set, for z G F, H\°' 3o (z) = h io > jo (x(z)) and H^ k) (z) = h io ' jo (y(z)) 
- the subscript 1 standing for the identity element of the group W = (£,??}. 

• Next, we define Hf 30 , Hf 30 on and H l v (h3 °, H l v odo on r](F) by 

Vz G £ (F) : H*? do (z) = H^' 30 (£(*)) , H*°' io (z) = -Hf > jo (z) - h^f + x V° (*) , 
Vz G r/ (F) : H*° (2) = H?' jo (r}(z)) , ff*°> io (z) = -H^ (z) - h^' 30 + x l °y jo (z) . 

• Then, we define H£f°, H^ jo on £ V (F) and F^' jo , H^ 30 on t/£(F) by 

Vz G £r? (F) : H^' (z) = fl*"*> (£(*)) , H^ j0 (z) = -H^ jo {z) - K§f + s V° {z) , 
Vz G (F) : H l v f° (z) = Hi°> jo ( V (z)) , (z) = -H^f (z) - h^f + s V° (*) . 

• At last, we define Hi°J and on £t/£(F) = rj^rj(F) by 

VzG^(F) : flgf (z)=H%*(t{z)), H£l°(z)=Hl^(r 1 (z)). 

Therefore we have, for each of the six domains w(F) of the decomposition (jlip . defined 
two functions Hw' J0 and Hw ,:>0 ■ Then, as said at the beginning of the proof, we set 
H io 'i°(z) = Hi£' jo (z) and H io >i°(z) = Ht jo (z) for all z G w(F) and all w G W. 

With this construction, (fT2|) and ()13ap are immediately obtained. To prove (|13b|) . we 
can use the fact that it is possible to express all the functions Hw' 30 , Hw' 00 in terms only 
of Hi ' 3 °,Hi ' jo ,hiy o ,x io yi° : we give e.g. the expressions of H^f and on £r?£(F) : 
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We therefore obtain (fl3b|) . since with (fT3al) we get #i°' io (£r)£(z)) +^ 0,J0 (£r/£(z)) +/4°o° = 
x l °yi° (£f]£,( z )) f° r 2 ^ l 7 ?^^); an d since = {1, ^, 77, 77^, ^77, ^77^}. □ 



Theorem 8. For any i,j,io,jo > Oj 
-[z - l/z ]/Q x 



Ql-0,30 _ 

hJ 2m\p 1<0 2 -4pi,ipi _i]V2 



Jexp(j6»)[0,oo] 2 „..,-„,, 



Z 



dz 



x{z) l y(z) 3: 

where 9 G [27r/3,7r] and where we have set fi x = zq + 1/zq — \z\ + 1/zi] = 4(^4 — l)(xi — 

1)/(X 4 - Xl) < 0. 



Proof. We have already noticed that the generating function of the Green functions G %0M 
defined in ©, is holomorphic in {(x,y) € C 2 : |x| < l,\y\ < 1}. As a consequence and 
by using ( 
double integrals : 



the Cauchy formulas allow us to write its coefficients G* ^ as the following 



1 



[2m] 



{\x\ = \y\=l} 



Q (x, y) xiyi 



dxdy. 



Then, by using the uniformization ([7]), the location of the cycles {|x| = 1} and {|y| = 1}, 
see Picture O the residue theorem at infinity and the Cauchy theorem, we obtain that : 



Ql>0,JO 



1 

2m 



H ioJo t z \ + H io,h f z \ + />,jo _ i j t \ 

K _L \J_ °'° x (z)dz, (15) 

cxp(te)[o,oo] [dyQ(x(z),y(z))]x(z) i y(zy 



9 being any angle in [27r/3,7r] - [27r/3,7r] because on the one hand, it is not possible to 
take 9 > tt, since exp(27r)[0, 00] is a singular curve for H l °' JO , and on the other hand, it 
is not allowed to have 9 < 2tt/3, since exp(2z7r/3)[0, 00] is a singular curve for H 10 ' 10 , see 
Proposition [6j 

Then, from ([7]) and from the fact that d y Q(x(z),y(z)) = d(x(z)) 1//2 , we remark that 
we have 

x' W /< W <«W.,, W ) = ^-lA.l/fa + iA.-h + iA.]) 

In this way and by using (|13b|) in (fT5|) . we get (JHJ). □ 



3 Proof of Theorem 1 : asymptotic of the Green functions 

Beginning of the proof of TheoremQ]. For any 9 € [27r/3,7r], the function x{z) l y{zy 
is, on exp(z0)[O, 00], larger than 1 in modulus, see Picture [3l Moreover, it goes to 1 when 
(and only when) z goes to or to 00. This is why it seems natural to decompose the 
contour exp(i9)[0, 00] into a part near 0, an other near 00 and the remaining part, and 
to think that the parts near and 00 will lead to the asymptotic of G 1 ^ , and that the 
remaining part will lead to a negligible contribution. In this way appears the question 
of finding the best possible contour in order to achieve this idea ; in other words, it is 
a matter of finding the value of 9 for which the calculation of the asymptotic of (|14p on 
exp(z6*)[0, 00] will be the easiest, among all the possibilities 9 6 [27r/3,7r]. 

For this, we are going to consider with details the function x(z) l y(zy , or, equivalently, 
the function Xj/i( z ) = ln(x(z)) + (j/i) ln(y(z)). Incidentally, this is why, from now on, we 
suppose that j/i S [0, M], for some M < 00. Indeed, the function is manifestly not 
appropriate to the values j/i going to 00, for such j/i, we will consider later the function 
{%/ j)Xj/i( z ) = ^ n (. x ( z )) + m (y( z ))- Nevertheless, M can be so large as wished, and, 
in what follows, we assume that some M > is fixed. 
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Now we set Xj/ii z ) = Y2p>o u p(j /i) 2 ? ■> this function is a priori well defined for z in a 
neighborhood of 0. Moreover, with 0, we obtain that vo(j/i) = and that for all p > 1, 

pv p (j/i) = {zl + 1/4 - [z{ + l/z[]) + (j/i)(4 + 1/4 - [zl + l/z p 3 ])/KP. (16) 

Likewise, we easily prove, by using ((TJ), that for z near oo, Xj/i( z ) = J2 p >o'^pU / ^) z P ■ 
Consider now the steepest descent path associated with Xj/h that is the function Zj/i(t) 
defined by Xj/i( z j/i(t)) = t- By inverting the latter equality, we immediately obtain that 
the half-line (l/i/i(j/i))[0, oo] is tangent at and at oo to the steepest descent path. 

Now we set, for the sake of briefness, pju = l/v\(j/i). With this notation, let us 
now answer the question asked above, that dealt with finding the value of 9 for which 
the asymptotic of G* ^ will be the most easily calculated : we are going to choose 
6 = aig(pj/ i ) € [27r/3,7r], and the decomposition of the contour exp(i#)[0, oo] will be 

exp(i0)[O,oo] = (p i/i /|^- /i |)([0,e]U]e,l/e[U[l/e,oo]). 



By using this decomposition in (|14|) . we consider now that the Green functions are the 
sum of three terms, and we are going to study successively the contribution of each of 
these terms. 

But first of all, we simplify the expression of pjn. Setting Sly = Z2 + I/22 — [ z 3 + V^] = 
4(y4 — l)(yi — 1) / (2/4 — y%) and using (fTBj) . we immediately obtain that v\(j/i) = 

+ {j/i)Qy/K- But it turns out that for all the walks of & a 8, we have Sl y = aSl x 
- this follows from a direct calculation starting from the explicit expression of the branch 
points xi, X4, yi, 2/4 in terms of the jump probabilities (pij)ij and by using Remark [H 
Therefore we have : 

111 

Pj/i ~ Mj/i) ~ (1 } 

Contribution of the neighborhood of 0. In order to evaluate the asymptotic of 
the integral (fl"4"|) on the contour (pj/i/\pj/i\)[0,e], we are going to use the expansion of the 
function (1/z) Ylw&wi~ l) 1 ^ x 1,0 y 30 (w(z)) at - expansion that we will obtain in Equation 
l\'21\) below. This is why we begin by studying the asymptotic of the following integral, for 
any non-negative integer k : 



z k 



-dz. (18) 



P S /i/\Pi/i\)[0,e] x(z) l y(z) 

By using the equality l/[x(z)' l y(z) 3 ] = exp(— iXj/i( z )) as wen as the expansion (fl~6]) of 
Xj/i at and then making the change of variable z = pj/%t, we obtain that (fl~8l) is equal to 

p)f f /lPj/l1 t k exp (-it) exp ( - iv 2 (j/i) ( Pj/i tf) exp ( - i £ u p (j/i) (Pj/itf) dt. (19) 

^° ^ p>3 ' 

Now we set m = max{|zo|> l/|-zo|> \ z i\i Vl^il) \ z 2\-> Vl^l, \ z z\i 1 / 1 ^3 1 } • Then, with (fl~6j) , 
we get \v v (j/%)\ < AmP(l + M). Thus, for all t € [0,e/\ Pj/i \], \ - i £ p > 3 ^(i/*)(Pj/i*) p l < 
4(1 + M)i(me) 3 /(1 - me). This is why exp(-i J2 P >3 t/ P U/ i ){Pj/i t ) p ) = 1 + 0(ie 3 ), the O 
being independent of j/i £ [0, M\ and of t G [0, e/\pj/i\}. The integral (fl"9|) can thus be 
calculated as 

(p j/t /i) k+1 [1 + 0(ie 3 )] r /lPj/ ' 1 t k exp (-t) [1 - v 2 (j/i) plf/i + 0(t*/i 2 )] dt. 

J 
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Let us now choose e such that ie/\pju\ — > oo and 0(ie 3 ) = o(l/i), e.g. e = 1/i 3 / 4 . For 
this choice of e, we obtain that the integral (|18|) is equal to 

(p 3/l /i) k+l [1 + o(l/i)] [k\ - v 2 (j/i) p) tl (k + 2)\/i + 0(1/* 2 )] , (20) 

where the o and O above are independent of j/i € [0, M\. 

We are presently ready to find the asymptotic of the integral (|14p on the contour 
(pj/i/\pj/i\)[0, e]. To begin with, we have the following expansion in the neighborhood of 
(directly obtained from (J7J), ([9]) and Remark HJ) : 

(-l) Hw) *V° (w(z)) = (z3 3 / 2 /2)a^ jo (fo + + /3) z 3 + O (z 6 ) . (21) 

Equation ([2T]) implies then that the integral (fl~4"|) on the contour (Pj/i/|Pj/i|)[0, e] equals 
-[zo-l/zo]/^ /• (^^a^oio^o + aio + ^^ + O^ 5 ) 



27rz[p lj0 2 - 4^1,^1 _!]V2 7( PjA /| PjA |)[ , e ] x (zy y (z) j 



-dz. 



So, with (|18p and (|2U|) applied for k = 2 and A; = 5, we obtain that the integral (|14|) 
on the contour (pj/i/\pj/i\)[0, e] is equal to 

-[z -lAo]3 3 / 2 a^ + + 3 _ 2 ^ 2 + 

Contribution of the neighborhood of oo. The part of the contour close to oo, 
(pj/i/\pj/i\)[l/e, oo], is related to the part (pj/i/\pj/i\)[0, e] by the transformation zt-^l/z. 
Moreover, it is immediate from ([7]) that for / = x, f = y, or / = Ylw<=w(~ ^y^x %0 y^°(w), 
f(l/z) = f{z). Therefore, the change of variable z H )• 1/z immediately entails that the 
contribution of the integral (jl4H near oo is the complex conjugate of its contribution near 0. 



Contribution of the intermediate part. We first recall from Proposition [6] that D 
denotes {xexp(i9) : x > 0, 2n/3 < 6 < ir}, and we define A e = {z € C : e < \z\ < 1/e}. 
Clearly (see Picture [3]) there exist r] Xj€ > and rj y ^ > such that for all z € D n A e , 
\x(z)\ > 1 + and |y(z)| > l + rjy te . In fact, since x'(0) = Q x Q and y'(0) = Q y /K ^ 0, 
it is possible to take rj Xt€ > rje and r\ y ^ > rje, for some 77 > independent of e small enough. 
Let us now consider the function 



s(z] 



x i( >yio(z) 



(-l) l ^x io y jo (w{z)) 



and let us show that sup zGD \s(z)\ is finite. For this, it is sufficient to prove that s has no 
pole in the closed domain D U {00}. 

By ©, the only zeros of the denominator of s are at z%, l/zi,Kzs,K/zs which, as we 
easily check, belong to —(DUD). Also, by (J7|) and ([9]), the only poles of the numerator of 
s are at K 2k z , K 2k / z , K 2k+l z 2 , K 2k+1 /z 2 , for k G {0, 1, 2}. Next, we verify that both z 
and Kz 2 belong to D, so that among the twelve previous points, in fact only z$ and Kz 2 
are in D. But in the definition of s, we took care of dividing by x^y^ , so that s is in fact 
holomorphic near these two points. Moreover, s is clearly holomorphic at 00. Finally, we 
have proved that the meromorphic function s has no pole in the closed domain D U {00}, 
hence s is bounded in D U {00}, in other words sup z£D \s(z)\ is finite. 

In particular, the modulus of the contribution of the integral (|14p on the intermediate 
part (Pj/i/\Pj/i\)]e, l/e[c D n A € can be bounded from above by 

\zp - l/zp\/\Q x \ 1 sup z£g |g(z)| 
2vr|pi i0 2 -4pi i ipi_i| 1 /2 e 2(i +r?e )i-io(i + 7?e )i-io- 1 > 
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Note that the presence of the term 1/e 2 in (j23|) is due to the following : one 1/e appears 
as an upper bound of the length of the contour, the other 1/e comes from an upper bound 
of the modulus of the term 1/z present in the integrand of (|14p . 

Then, as before, we take e = l/i 3//4 , and we use the following straightforward upper 
bound, valid for i large enoug h : 1/(1 + r ? /i 3 / 4 ) i < exp(-(r ? /2)i 1 / 4 ). We finally obtain 
that for i large enough, ([23]) is equal to 0(i 3//2 exp(— (ry/2)f 1//4 )). We are going to see soon 
that this contribution is negligible w.r.t. the sum of the contributions of the integral (|14[) 
in the neighborhoods of and oo. 

Conclusion. We have shown that the contribution of the integral (|14p in the neigh- 
borhood of is given by (I22p . that the contribution of (I14p in the neighborhood of oo is 
equal to the complex conjugate of (|22p . and that the contribution of the remaining part 
equals 0(i 3 / 2 exp(— (ry/2)i 1 / 4 )). Moreover, starting from (117p . we immediately get that 
(Pj/i/i) 3 - (~pjji/i) 3 = -i3 3 / 2 aij(i + aj)/[U x (i 2 + aij + a 2 j 2 )] 3 . In this way, we obtain 



in -\z - l/z ]3 3 / 2 aQl .... . n . 
G i = A T 2 / rn^oJo^o + aj + P x (24) 

^3 3 / 2 mj(* + aj) ^U/^Pj/i 5 ~ ^U/i)pjJi 5 { ^/j^ 



[Vt x {i 2 + aij + a 2 j 2 )Y 



If 7 G]0,7r/2[ and j/i — > tan(7), then + aj)/(i 2 + aij + a 2 j 2 ) 3 ~ C 7)Q ,/i 3 with 
C 7jQ > : Theorem [T] for 7 g]0, tt/2[ is thus an immediate consequence of ([2"4"]) . In that 
case, there was in fact no need to make an expansion with two terms in (|22p and (|24l) 
above, one single term would have been accurate enough. 

If j/i — > tan(O) = 0, then + aj)/(i 2 + aij + a 2 j 2 ) 3 ~ (j/i)/i 3 . By using the 
explicit expressions of ^2 OA) an d /Oj/i, see respectively (fTBT) and (fTTj) . we easily obtain 
that z^2(j /i)pj/i 5 — ~&2(j / i)j>JJi 5 = 0(j/i). This implies that the sum of the two last terms 
in the brackets of (f24l) equals 0((j /i)/i 4 ) + o(l/i 4 ), which is obviously negligible w.r.t. 
(j/i)/i 3 . Theorem [1] is therefore also proved in the case 7 = 0. 

In order to prove Theorem [1] in the case 7 = vr/2, we would consider (i/j)Kju rather 
than Kj/i, and we would use then exactly the same analysis, we omit the details. 

Finally, in order to prove that the constant C in the statement of Theorem[T]is positive, 
it is clearly enough to show that i[zq — l/zo]/[(pi,o 2 — ^Px^Pi-i) 1 ^^] is positive. 

For this, note first that from its definition, it is immediate that VL X < 0. Moreover, 
it follows from the beginning of Section [2] that if X4 > 0, then i[zq — 1/zq] < and 
(Pi,0 2 -4pi,ipi,_i) 1 / a > ; if x 4 < 0, then [z -l/z ] > and i/(pi, 2 -*Pi,\Vi ,-i) 1/a < ; 
and if X4 = 00, by taking the limit in anyone of the two previous cases, we obtain that 
i[z - 1M)]/[Cpi, 2 - 4p ljl p lj _ 1 ) 1 / 2 ] < 0. 



A few words about the analytical approach used here. The two key steps in the 
proof of Theorem Q] are first the explicit expression for the Green functions (I15p , and then 
the expansion (il3b) of H io ' jo + H io > jo + h^J° - x io y jo at 0, which is the numerator of the 
integrand in f)15|) . 

It is worth noting that for any walk of & D it is still possible to obtain (|15f> 

- without additional technical details, besides. On the other hand, obtaining explicitly 
the expansion at of H l0 > 30 + H l0 > 30 + Hqq — x^y 30 in the general setting seems us quite 
difficult - all the more so as this expansion has to comprise several terms, since a priori 
it could happen that several terms lead to non-negligible contributions in the asymptotic 
of the Green functions. 

It is more imaginable (though technically difficult) to obtain this expansion for the 
walks for which an equality like (|13bp holds ; unfortunately, having such an equality is far 
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for being systematic, even for the processes associated with a finite group W : for example, 
the random walk with jump probabilities pi t i = Po,-l = P-1,0 = 1/3 has manifestly a group 
W of order six, but does not verify an identity like (|13b|) . 
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